Impact of observational uncertainties on universal scaling of MHD turbulence 
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Scaling exponents are the central quantitative prediction of theories of turbulence and in-situ 
satellite observations of the high Reynolds number solar wind flow have provided an extensive 
testbed of these. We propose a general, instrument independent method to estimate the uncertainty 
of velocity field fluctuations. We obtain the systematic shift that this uncertainty introduces into 
the observed spectral exponent. This shift is essential for the correct interpretation of observed 
scaling exponents. It is sufficient to explain the contradiction between spectral features of the 
Elsasser fields observed in the solar wind with both theoretical models and numerical simulations 
of Magnetohydrodynamic turbulence. 
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Universality in isotropic, homogeneous turbulence is 
expressed through its statistical scaling properties. In the 
absence of intermittency, the scaling exponent for the in- 
ertial interval of hydrodynamic turbulence is completely 
determined by the assumption of self-similarity [H , lead- 
ing to the well known unique —5/3 power spectral ex- 
ponent. This is not the case for Magnetohydrodynamic 
(MHD) turbulence where the magnetic field introduces 
an additional physical quantity with the dimension of 
velocity, namely, the Alfven velocity Q and indeed it is 
an open question as to whether the scaling is universal. 
Detailed phenomenological models of MHD turbulence 
are thus needed to predict the scaling exponent, and its 
precise observational determination is essential in order 
to validate these theories. 

In-situ satellite observations of the solar wind mag- 
netic field and bulk flow span several decades in tempo- 
ral scales and offer a 'natural laboratory' for the study 
of MHD turbulence, the Reynolds number exceeds ~ 10^ 
Q. They have been extensively used to test theoretical 
predictions of MHD turbulence (see, and references 
therein). The Elsasser fields, = v± 'B/y/Anp, where 
V and B are the velocity and magnetic fields, respec- 
tively, and p is the average density, represent eigenfunc- 
tions of counter propagating (with respect to the mean 
magnetic field) Alfven waves and therefore they are pri- 
mary fields for the study of incompressible MHD turbu- 
lence. Fluctuations in the fast solar wind are strongly 
imbalanced - there is more power in Alfven waves prop- 
agating outward from the sun than toward it (e.g., |ig |) 
so that the power in Z+ dominates over that in Z~. 
As with many other quantities that characterize phys- 
ical properties of the turbulent flow (e.g. Yaglom re- 
lations |lll4l3t and dynamic alignment angle [7|), the 
Elsasser variables combine velocity and magnetic field 
fluctuations as a function of temporal scale. Pioneering 
observations from the HELIOS missions showed that in 



the fast solar wind streams at 1 AU the observed power 
spectrum of Z~ (the subdominant component) did not 
follow a single power law shape. At very low frequen- 
cies (/ < 3 X 10~*Hz) the spectral slope 7_ was close 
to Kolomogorov's value (7_ « —1.67), whereas at higher 
frequencies (5 x 10~*Hz < / < 2 x 10^'^Hz) the Z~ power 

fectrum was much shallower, with 7_ w —1.3 — 1.4 0- 
A similar trend was found more recently in WIND 
observations Q where in the low frequency part of the in- 
ertial interval (10~^Hz < / < lO^^Hz) Z~ nearly follows 
Kolmogorov scaling which at higher frequencies again is 
more shallow (7_ « —1.3). The absence of single scaling 
of the subdominant Elsasser field in the inertial inter- 
val contradicts all recently developed models of strong, 
anisotropic imbalanced MHD turbulence 14-1^ which 
predict a universal scaling for both dominant Z"*" as well 
as sub-dominant Z~ spectra. They also are inconsistent 
with the results of recent high resolution direct numer- 
ical simulations of imbalanced MHD turbulence which 
showed nearly the same spectral indices of the energy 
spectra in the inertial interval 15l. Il7|. 



Control of observational uncertainty in the in-situ ob- 
servations is non-trivial, although these errors often have 
known bounds. There are different challenges for mag- 
netic field and velocity measurements; solar wind velocity 
observations are intrinsically more uncertain compared to 
the magnetic field data [19|. In this Letter we propose 
a general, instrument independent method to estimate 
the uncertainty on velocity field fiuctuations direct from 
the data. We obtain the systematic shift that this uncer- 
tainty introduces into observed spectral exponents. We 
will see that the shallower Z~ spectrum at high frequen- 
cies can be entirely accounted for by this uncertainty in 
the velocity data and the observations of the Z^ spectra 
may in fact within achievable accuracy of the observa- 
tions, be in agreement with the predictions of theory and 
numerical simulations. 



2 



We use data obtained by the WIND spacecraft at 3 
second resolution. Magnetic field data is provided by 
the MFI instrument [20|jand density and velocity data 
by the 3DP instrument [2l|. We use observations made 
during a quiet fast stream of April 04-06, 2008, during 
which the solar wind speed remained above 550 km/s. 
The energy of compressive fluctuations was an order 
of magnitude lower than that of incompressible fluctu- 
ations and, consequently, magnetic and velocity fluctu- 
ations, being mainly Alfvenic, were dominated by the 
components perpendicular to the local mean field. The 
mean field, B(i, r), at some time t and on scale r is 
defined as the magnetic field averaged over the interval 
[t — T,t + 2t] . Similar to most studies of the Alfvenic 
component of fluctuations (e.g., [zl) we focus on the per- 
pendicular components of the fluctuations of the velocity 
defined as 6v± = (5v — (B • Sv)Sv and magnetic field 
SBj_ = SB-(B SB)SB where Sv{t, r) = v(t + t) - v{t) 
and SB{t,T) = B(t + T) — B(i); in what follows subscript 
_L will be omitted for simplicity. 

In common with all velocity in situ observations, the 
3 s velocity observations on WIND (as well as any other 
in-situ measurements of moments of the particle distribu- 
tion function) are quantized before ground transmission 
and this quantization results in high frequency noise or 
quantization noise jlQi] . This, and other contributions to 
observational uncertainty, decorrelate the velocity and 
magnetic field fiuctuations at high frequencies. White, 
delta correlated noise provides a reasonable generic, in- 
strument independent model for the uncertainty [1] . Any 
measurement of a velocity component fiuctuation Svo can 
then be represented as a sum of the 'real' turbulent signal 
Svg and a noise 5vn which has zero mean and standard 
deviation e, so Svo = Svs + Svn- In what follows we will 
neglect the uncertainties in the magnetic field measure- 
ments since generally these are small relative to that on 
the velocity measurements [1, |^ . 

We will first quantify the velocity uncertainty from the 
data. We will exploit the fact that both the turbulent 
signal and the noise are random variables with distinct 
characteristic autocorrelation time scales. We make a key 
assumption- that the autocorrelation timescale of the un- 
derlying turbulent signal is that observed in the magnetic 
field component fluctuations SBo (they have negligible 
noise) and that this is also the autocorrelation timescale 
of the 'true' turbulent velocity component fluctuations 
6vs. Any difference in the autocorrelation functions of 
the observed 6vo and SBo are thus attributable to the 
(delta correlated) noise on the velocity Svn . The autocor- 
relation coefficient (AC) of a component 5v on time lag A 
is RsAr^A) = E[iSv{t+A,T)-{Sv))iSv{t,T)-{Sv))]/al, 
where ] is the expected value operator and asv is the 
standard deviation on a given velocity component. The 
autocorrelation coefficients i?5B^(r, A) and i?5„^(r, A) 
are plotted for lag A = 3s as a function of scale r with 
black and red lines respectively in Figure 1. We see that 
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Figure 1: Autocorrelation functions RsBy{T, At) (black line), 
Rsvy {t, At) (red line) and Rsvp {t, At) (dashed line), with the 
time lag At = 3 s. 



the AC grows with scale r for both signals and that the 
velocity AC is systematically lower that that of the mag- 
netic field, consistent with the assumption of delta or un- 
correlated noise (< 5vn{t + A)Svn{t) >= (5(A)) that prin- 
cipally affects the velocity signal. Given these assump- 
tions one can construct a pseudo noisy signal by adding 
uncorrelated noise to the magnetic field observations. 
The pseudo noisy signal fluctuations 6Bo+n ~ SBo+SBn, 
where SB^ are delta correlated Gaussian distributed ran- 
dom numbers with zero mean and standard deviation 
Eb- The magnitude of the pseudo noise eb can then 
be systematically varied and we plot on Figure 1 (black 
crosses) the result for a fractional uncertainty on the 
magnetic field corresponding to a velocity uncertainty of 

£b = £b^ (^"o)/('5-So+ri) = 4 km/s. We see that this 
pseudo noisy signal closely coincides with the observed 
velocity AC, suggesting that e ^ 4 km/s is a reasonable 
estimate of the amplitude of the noise on the turbulent 
velocity signal. We will develop this idea to obtain gen- 
eral methods to estimate the uncertainty direct from the 
data. First, we will see how these uncertainties can affect 
measurements of scaling exponents and the conclusions 
that can be drawn from them. 

In Figure 2 we plot the observed second order struc- 
ture functions 5*2 of a component of the fluctuations in 
the Elsasser variables, where the Elsasser components are 
given by (5Zj^^(r) = 5vi{t, T)±SBi{t, r)/-\/47rp and p(r) is 
the interval averaged local mean value of the density over 
time scale r, and = {6Zf{T)'^). The solid lines are 
the structure functions of z components of the dominant 
(black) and subdominant SZ~ (red) fields. They 
are normalized to have the same values at r = 10 min 
scale on this plot; the power in 6Z^ is 20 times that in 
SZ~ . For the ideal statistical scaling of fully developed 
MHD turbulence we anticipate the scaling 5*2 ^ r'»±'^^^ 
and turbulence theories predict constant values of C^(2) 
over the entire inertial interval (they are directly related 
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Figure 2: The normalized second order structure functions: 
of the 2 component of the sub-dominant Elsasser variable 
(red solid line), of with added Gaussian noise with 
£fl = 4 km/s (red crosses), of the dominant Elsasser variable 
Z'^ (blue solid line) and Sa (circles, see text for details). In 
the insert: the second order structure functions of the sub- 
dominant Elsasser variable for different values of added Gaus- 
sian noise. Raw observations are denoted by the solid line, and 
with added noise e_g = 2 km/s (dashed line), es = 4 km/s 
(dash-dotted line) and eb = 5 km/s (dotted line). 



to the power spectral exponents 7± via 7± = — C±(2) — I)- 
We can see that, consistent with earlier studies [^-Isllsl, 
the subdominant Elsasser variable does not follow a sin- 
gle power law in the inertial interval. A linear mean 
least square fit over scales 30 s < r < 10 min gives 
7+ = -1.54 ± 0.02 and 7_ = -1.40 ± 0.02, consistent 
with previous observations. 

We now test the idea that the uncertainty in the ob- 
served velocity, estimated above, is sufficient to account 
for this observed departure from the theoretical predic- 
tions. Since the power in JZ"*" is significantly higher than 
that in 6Z^ we will focus on the effects of uncertain- 
ties in velocity on the SZ~ signal only. We calculate 
S2 iSBo+n,Svo) using the pseudo noisy signal fluctua- 
tions SBo+n and the observed Svq. This is shown in the 
inset of Figure 2 for a range of values of £_b. We can see 
that addition of 'white' (delta correlated) noise always 
systematically 'flattens' these curves, that is, it decreases 
the value of the scaling exponent; for — A km/s pseudo 
noise strongly affects S2 at all scales in the inertial inter- 
val. The 'flattening' of the pseudo-noisy S2 {SBo+n,^Vo) 
curve, that is, the change in the mean exponent over 
timescales 30 s < r < 10 min is A7_ « 0.13, is close to 
the observed difference between exponents of the domi- 
nant and subdoninant fields (7- — 7+ = 0.14) hence this 
difference could be just due to noise in the velocity data. 

The pseudo noisy S*^ curve generated with noise es = 
4 km/s is plotted as the red crosses in the main panel 
of the Figure. Since the noise is assumed to be linearly 
additive, the difference between the observed, and the 
pseudo noisy S2 curves, that is, = [S2 (SBo+m^Vo) — 
S2 {SBo, Svo)] provides an estimate of how a velocity un- 



certainty of £ = 4 km/s affects the subdominant El- 
sasser variable scale by scale. We then compensate 
for this systematic effect by subtracting this r depen- 
dent uncertainty from the observed S2 [SBq, 6vo) and 
S2 {5Bo,5vo) — es is shown by the black circles on the 
plot. This compensated subdominant 5*2^ curve now has 
a single scaling range, consistent with current theories 
and numerical predictions 14l4l8l| . It can also be seen to 
be in remarkable agreement with that observed for the 
dominant Elsasser variable. The uncertainty in the ve- 
locity that we have estimated from the data, as shown in 
Figure 1, is thus sufficient to account for the departure 
in scaling between the 5Z^ and (5Z+ Elsasser variables 
and these observations may in fact within the achievable 
accuracy be in agreement with theories [H, 17, 1^ that 
predict a single scaling for SZ~ and JZ"*". 

We have verified that (5Z+ is not strongly affected by 
velocity uncertainties of this amplitude: S2 remains al- 
most unchanged for sb < 5km/s. This conclusion is also 
supported by the fact that as shown in Figure 2, S2 
has a convex shape, which is typical for finite range hy- 
drodynamic and MHD turbulence 13, 23 1, whereas S2 
is concave, consistent with 'shallowing' at small scales 
due to noise effects. In addition to the velocity measure- 
ment uncertainties the Elsasser fields are also affected by 
uncertainties related to the density measurement. Our 
analysis shows that adding the same amount of relative 
pseudo noise to the density data had negligible infiuence 
on the results. This is not surprising, as the definition 
of the Elsasser fields is in terms at the mean value of 
the density over scales, and this reduces the influence of 
density uncertainties due to the central limit theorem. 

We now outline methods to explicitly determine the 
uncertainty as a function of scale t from the observations. 
We will generalize the approach shown in Figure 1. Given 
the assumption that i?5B^(T, A) = Rsv^{t, A) and again 
that the noise is delta correlated < Svn{t+A)Svn{t) >= 
we can obtain the uncertainty from the AC directly: 



i?fa„(T,A) 



(1) 



Alternatively, we can estimate the uncertainty in velocity 
fluctuations using a pseudo noisy signal constructed by 
adding uncorrelated noise to the magnetic field observa- 
tions: 5Bo+n — SBo + SBn- We can vary Sb until the 
AC of the pseudo noisy magnetic field RsBa+r, ("^j ^) co- 
incides with that of the velocity Rsy^ (r. A). The relation 

eB{A.)/ ^ {5Bl^^) ~ £{A)/ ^ {Sv'^) provides a scale de- 
pendent estimate of e. These methods are demonstrated 
in Figure 3 and give results that are consistent with the 
estimate of e = 4 km/s used above. 

Recently, in-situ solar wind data have been used 0] to 
test the prediction of scale dependent dynamic alignment 
in MHD turbulence 2J|. This relies on determination 



of the angle between magnetic and velocity fiuctuations 
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perpendicular to the local mean magnetic field direction. 
The observational uncertainty is known to make a sig- 
nificant contribution to the component of the velocity 
perturbation perpendicular to the magnetic field pertur- 
bation S\'p{t,T) = (5Bj_ X Sv±/SB± even at quite large 
scales T ^ 10 min We show this in Figure 1, where 
the AC of Svp{t,T) is plotted (dashed line), we can see 
that the magnitude of the AC is much lower than that 
of the components of velocity and magnetic field. This 
AC function allows us to develop one more method to 
determine measurement uncertainty. Representing Svp,o 
as a sum of the underlying turbulent signal and a noise 



Svr, 



Svp,n and assuming that the AC of Svp,. 



and the magnetic field fluctuations are identical, we have 
that: 



(2) 
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Figure 3: Quantization error derived by different methods (see 
text for details): green line correspond to the result derived 
usimg Eq. U]), blue line is derived by adding artificial noise 
to the magnetic field data. Red line corresponds to the result 
derived using Eq. ([2]). Black line denotes rms value of Svp. 



The scale dependent estimate of e derived using this ap- 
proach is given by the red dash-dotted line in Figure 3 
and the result is in very good agreement with the results 
obtained by the other methods. We also plot the rms 
value of the observed Svp^o on Figure 3 (black line) and 
this can be seen to coincide with our various estimates 
of the uncertainty in the velocity fluctuations on scales 
T of a few minutes or less. On these smallest scales, the 
observed (5vp_o is almost entirely dominated by observa- 
tional uncertainty. 

The estimate of the quantization error in [tI gave 
a somewhat lower value e « 1.5 km/s. In [7[ the 
alignment angle 9a = aicsm{{6vp6B±) / {Sv±6B±)) « 
Vpl \J {6v±{t)^) was used to estimate the quantization 
error. In fact correlation between SBj^ and 6v±^ is much 
stronger than correlation between 6B± and Svp (because 
the latter at small scales is strongly dominated by the 
error). Detailed analysis shows that this increases the es- 
timates in Q by a factor ~ 1.6. Also, the local alignment 
angle between velocity and magnetic fields 6 is random 
variable with zero mean. Any measure of the alignment 
angle (such as mean of the absolute value of the align- 
ment angle) should thus be related to the standard devia- 
tion ae of the angle. If for highly aligned cases we assume 
that 9 is close to Gaussian distributed, then its absolute 
value is described by the half-normal distribution and its 
mathematical expectation is given by ag y'L — 2/Tr (and 
not ae) and this yields another factor of 1.2533 compared 
to the estimate of {3] ■ 

This uncertainty of « 4 km/s in the velocity fluctua- 
tions will systematically reduce the scaling exponent of 
the velocity compared to that of the magnetic field, if as 
we have assumed here, the uncertainty in the magnetic 
field can be neglected in such a comparison. We estimate 
that this uncertainty alone would result in a difference in 
the power spectral exponents of — 1b ~ 0.04. This 
is less than the observed difference which is typically in 



the range ~ 0.1 [17[ to ~ 0.2 [25|- These estimates are 



sufficiently close however to merit more detailed analysis. 

In summary, we have presented general, instrument 
independent methods to determine the uncertainty in 
velocity fluctuations in single point measurements. We 
have shown that this uncertainty is sufficient to account 
for the departure in scaling between the subdominant 
and dominant Elsasser variables and thus are able to re- 
port for the first time that the observations are, within 
the achievable accuracy, in agreement with the predic- 
tions of theory and numerical simulations. Our results 
put careful estimation of uncertainties at the centre of 
the testability of theoretical predictions of scaling expo- 
nents. Our approach, and development of it, is highly 
relevant for methods developed for the study of MHD 
energy cascade rates in the solar wind. These inevitably 
involve combinations of velocity and magnetic field fluc- 
tuations that are scale dependent such as_ mixed third 
order moments of the Elsasser fields 
relations flU. 
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